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Abstract 

A unified QCD formulation of leptoproduction of massive quarks in charged current 
and neutral current processes is described. This involves adopting consistent factor- 
ization and renormalization schemes which encompass both vector-boson-gluon-fusion 
( "flavor creation" ) and vector-boson-massive-quark-scattering ( "flavor excitation" ) pro- 
duction mechanisms. It provides a framework which is valid from the threshold for 
producing the massive quark (where gluon-fusion is dominant) to the very high en- 
ergy regime when the typical energy scale is much larger than the quark mass mg 
(where the quark-scattering should be prevalent). This approach effectively resums all 
large logarithms of the type (Q;s(/i) log(/i^/mQ))" which limit the validity of existing 
fixed-order calculations to the region /i ^ 0{niQ). We show that the (massive) quark- 
scattering contribution (after subtraction of overlaps) is important in most parts of the 
{x, Q) plane except near the threshold region. We demonstrate that the factorization 
scale dependence of the structure functions calculated in this approach is substantially 
less than those obtained in the fixed-order calculations, as one would expect from a 
more consistent formulation. 
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1 Introduction 



The production of heavy quarks in photo-, lepto-, and hadro-production processes has be- 
come an increasingly important subject of study from both the theoretical and experimental 
points of view. However, there are some outstanding problems with existing perturbative 
QCD calculations of heavy quark production: sizable (spurious) scale dependence of the 
predictions, apparent disagreement with observed b-production cross-section at the Teva- 
tron, ... etc. See Ref.|l| for a recent review of the theory and phenomenology of heavy 
quark production. As will be discussed later in this paper, there is also an inconsistency in 
most theoretical calculations of the cross sections: the schemes used in the next-to- leading 
order calculations are not the same as those used in determining the parton densities from 
global analyses, such as in Refs. nil. In this paper we will spell out the details of 
a more complete and consistent formulation of heavy quark production. For the sake of 
clarity, we shall focus on the case of leptoproduction, although the same principles apply 
to hadroproduction as well. 

Perturbative QCD calculations rely on factorization theoremsEC: Different factors in- 
volve different scales of virtuality, and a factor that involves only physics on a scale m 
can be effectively calculated in a power series in as{m). The simplest factorizations, like 
the operator product expansion, are for certain two-scale problems: One scale, about the 
bound-state nature of hadrons, is of the order of A or the mass of a typical hadron; and 
the other, defining a scale of large virtuality, can be in deep-inelastic scattering or the 
transverse momentum of a measured jet]^ 

Processes involving heavy quarks are a good example of a multi-scale process, for in 
practice we may then have to deal with at least four scales (which we denote by fii): A, 

(as above), and the masses of the charm and bottom quarks, rric and mi,. When one 
uses conventional calculational schemes designed for two-scale problems, the presence of 
more than one large scale results in logarithms ln(/ij//ij) of the large ratios in the higher 
order correction terms. These logarithms vitiate the very basis of the original perturbative 
calculation, because of the large size of the yet higher order terms beyond the order included. 
In view of the high energies at present colliders, this problem can defeat the large effort put 
into existing fixed order calculations. 

The basic principles for treating this situation were constructed a long time ago by 
Witten in his workM on heavy quarks in the operator product expansion (OPE). But his 
methods as they stand do not provide a sufficient algorithm for calculating the processes in 
which we are interested. An important practical concern is that the methods of calculation 
should be applicable when some of the scales of interest are comparable to each other as 
well as when they are very different. Simple minded methods involYing "integrating out 
the heavy quarks" are not sufficient, but, as Collins, Wilczek and ZeelH pointed out a long 
time ago, the problem can be conveniently considered as one of choice of the subtraction 
scheme (for renormalization and factorization). We shall refer to this work as CWZ in the 
following. 

Consider the case of deep-inelastic scattering, with Q denoting the invariant mass of the 

^To avoid circumlocutions, we will often use the terminology of the operator product expansion when 
discussing factorization theorems. In particular, we will use the term "Wilson coefficient" to denote the 
short-distance coefficient in the standard factorization theorem for deep inelastic scattering etc. 
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exchanged boson, and with one heavy quark, of mass m. Three kinematic regions are of 
interest: 



• m ^ Q: The quarkmass is larger than ah other scales in the problem, so that the 
decoupling theoremO applies: All graphs involving the heavy quark may be dropped, 
at the price only of a possible finite renormalization of the parameters of the theory, 
notably as- 

• m = 0{Q): The heavy quark mass must be treated in the same way as Q — as 
a large parameter. Heavy quark lines appear in the Wilson coefficients or in finite 
renormalizations of exactly the same kind as in the decoupling theorem. 

• m ^ Q: As far as the OPE is concerned, the heavy quark is to be treated as light: its 
mass is to be neglected in Wilson coefficients, and there are parton densities for the 
heavy quark. Since the quark is heavy on an absolute scale, i.e., m » A, Witten's 
methods may be used to calculate its density in terms of the densities of light partons. 

The method we will describe will give a unified treatment that will cover all ranges of 
mass. Furthermore, the CWZ method also allows us to treat the case that there are several 
heavy quarks, whose masses may or may not be strongly ordered. 





Qi 



Qi 3^ 



Q2 



Q2 



Qi 



(b) 



Figure 1: Amplitudes for heavy quark production: (a) order quark-scattering; and (b) 
order a\ gluon- fusion contributions. At least one of the quarks, say Q2, is "heavy" and 
corresponds to the "Q" used in the text. For fiavor changing currents, the two quarks Qi 
and Q2 are different. For neutral currents, they are the same. 

One can see the issues in contrasting treatments of heavy quark production in lepton- 
hadron scattering in existing literature. For charged current interactions, such as charm 



production in neutrino processes, most existing work focuses on the dominant underlying 
order-Qfg parton pro cess W + s ^-Lpf- Fig- the quark-scattering or "ffavor-excitation" 
sub-process).'— ' HI' — ' — ' 0' — ' In contrast, for neutral current interactions, such as 
charm and bottom production in electron- and muon-hadron scattering, practically all cal- 
culations begin with the order-a^ parton process 7/Z° + g ^ c + c {cf. Fig. the 
gluon- fusion or "ffavor-creation" sub-process). HI' HH In both cases the initiating parton is 
a light parton, and an appropriate scheme for computing radiative corrections is what we 
will call a 3-flavor scheme, where the heavy quark only occurs in the Wilson coefficients, 
and where there is no parton density for the charm and bottom quarks. 

But at very high energies (such as are now available at the HERA e — p collider and 
beyond) charm and bottom quark masses can become small compared to a typical energy 
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scale: we can have ^ m^, m?. It is then natural to count these heavy quarks, especially 



the charm quark, as partons-cOl c3 t3 We then effectively use a 4- flavor or 5-flavor 
scheme (as is done in commonly used parton distributions 0' B 0). In particular, the 
lowest order Wilson coefficient for charm production in the neutral current process is not 
gluon fusion, but flavor excitation + c — > c. Of course, the gluon density is rather 

large compared to the charm quark density, so that higher order gluon fusion process can 
be numerically comparable or larger than the flavor excitation process. (It is misleading to 
argue that an order al subprocess is smaller than an order subprocess merely by virtue 
of its being higher order unless the initiating partons are the same.) Moreover, one must 
make the correct subtraction from the gluon fusion process to avoid double counting. 

Evidently, for this purpose, the notion of a quark with mass mg being "heavy" or "light" 
must be taken as relative — with respect to the energy scale of the probe Uphy The latter 
view forms the very basis of the QCD parton model for the well-known light quarks u, d and 
s (which do have non-zero, albeit small, masses!). The parton approach effectively resums 
large logarithms in fixed-order calculations arising from initial state collinear singularities 
of the form (a^ log /ip/iy/mg)" to all orders in ag- 

On the other hand, one also has an absolute notion for the term "heavy quark" — that its 
mass is sufficiently large compared to A so that as{mQ), the effective coupling at the heavy 
quark mass, is in the perturbative region. This notion then refers to the charm, bottom 
and top quarks as heavy, regardless of the magnitude of the typical energy scale (Uphy) of 
the problem. This view is taken in all rie yt-t,o-leading order calculations on heavy quark 
production in the existing literature. lH' cB' c3 

As described earlier, the theoretical basis for a unified QCD treatment of heavy quark 
in the relative sense, suitable for all energy scales, already exists in the literature: it is based 
on the CWZ renormalization scheme which n atur ally implements the intuitive energy-scale- 
dependent light and heavy quark concepts. c2l This scheme has been applied with good 
success to the calculation of Higgs boson via heavy quarks and gluons, effectively unifying 
the corresponding quark-scattering and gluon-fusion subprocesses in one consistent scheme 
which is valid at all energies. This approach clearly also provides a natural framework 
for calculating the production of heavy quarks. It is particularly simple to implement in 
lepto-production processes, as already pointed out in a previous short communication. H^l 

The present paper presents details of the method and the main physics results. For defi- 
niteness, we shall refer to this approach as the variable (i.e., scale- dependent) flavor number 
scheme, in contrast to theHaxd flavor number scheme used in conventional calculations of 
heavy quark production. c3' c3 Since charged current and neutral current processes are 
treated in one uniform framework, we shall use the most general couplings for the vector 
gauge boson to the leptons and the quarks; and we shall keep the most general mass config- 
urations for the quark lines in the hard cross-section calculation. The resulting complexity 
in kinematics, in the application of the factorization theorem of QCD, and in the calcula- 
tion of hard matrix elements can be effectively handled using the helicity formalism. This 
aspect of the problem is formulated and presented in a separate paper, 1^ hereafter referred 
to as I. Sec. |2| provides an overview of the scale-dependent parton flavor number scheme. 
Sec. ^ gives the detailed results on the order gluon-fusion amplitudes. Sec. |^ discusses 
the subtraction procedure needed to make quark-scattering and gluon fusion mechanisms 
consistently co-exist. Sec. ^ presents the main physics results to show the relative impor- 
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Figure 2: Regions of the x — jjLphy kinematic plane for a typical physical process involving 
a heavy quark with mass mq and the natural QCD calculational schemes for each region. 

tance and interplay of the various contributions, as well as the reduced scale-dependence of 
the predictions. Finally, Sec. ^ recapitulates the theoretical issues and points to potential 
applications. 

Since our approach effectively resums all large logarithms of the type (a^ log fiphy/mg)"^ 
which occur in fixed-parton-flavor-number calculations, it naturally extends the range of 
validity of the latter beyond the region Hp^y ~ 0{mQ). It does not, however, deal with 
the class of logarithms of the type (a<j logs///^^^)" which is associated with the "small- 
x" problem (typically, x = fiphy/\/s)- The latter has been the subject of several recent 
studies and it r equ ires an entirely different method of resummation - the so-called kt- 
factorization.'H' These two approaches are compatible and complementary: they 
both extend the region of applicability of the perturbative QCD calculations, but to different 
regions of phase space.^ This is illustrated schematically in a map of the x — fi kinematic 
plane. Fig. |2[ Broadly speaking, our approach is needed when the typical energy scale 
fJ'phy becomes large (compared to mq) for not-too-small x; and the /ct-factorization method 
is necessary for very small-x and moderate fJ,phy What values of x must be considered 
as "small", and fiphy/mg as "large", to require these improvements are open questions 
with no easy theoretical answers in perturbative QCD. (c/., the similar question: for what 
value of Q should Bjorken scaling set in?) However, these questions can be investigated 
phenomenologically by comparing numerical results from the different approaches in their 
regions of overlap. Existing numerical studies of the small-x resummation and conventional 
approaches suggest the latter maybe valid down to x ~ 10~^.llll The results presented later 
in this paper will shed some light on the comparison of scale-dependent and fixed parton 
flavor number schemes. 



unified treatment is a topic for the future. 
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2 Overview of the Scheme and the Calculation 



We consider a general lepton-hadron scattering process: 

h{h) + N{P) ^2(^2) + Q{pq) + X{Px), (1) 

to lowest order in the electro-weak interactions, as depicted in Fig. ^. In the final-state, we 
have required that there be a heavy quark Q of momentum Pg-Q We label the exchanged 
vector boson (7, W, or Z) by B and its momentum by q. 




Figure 3: General lepton-hadron production amplitude for a heavy quark 



2.1 Hadron Structure Functions and Factorization 

After the calculable leptonic part of the cross section has been factored out, we work with 
the hadronic process induced by the virtual vector boson B: 

B{q)+N{P)^Q{pQ)+X{Px), (2) 
and the cross-section is expressed in terms of the hadronic tensor 

^"' = i^ E (PI-J'^IPq^ Px)i27r)^6^^^ (p + q-p^- Px){Px , PQ\r^\P) ■ (3) 

X{Px),spin 

where J2 denotes a sum over all hadronic stateq containing the final-state quark Q of 
momentum Pq. 

The factorization theorem asserts that the hadronic tensor has the form@ 

to the leading power of q^. Here /^(C,/^) is the distribution function of parton a in the 
hadron A'', and ka, as(/u)) is the Wilson coefficient. That is, Q'^^ is the same kind 

of object as the hadronic tensor Eq. (^) except that it is evaluated on a parton target, 

*We mean here a heavy quark in the absolute sense: mg 2> A, so that as{mQ) is smaU enough to be in 
the perturbative regime. The heavy quark will be detected by its hadronization products. 
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Figure 4: Graphical representation of the factorization formula, Eq. (|5|). 



and that the long-distance contributions are subtracted off. The proof of the factorization 
theoqrem is to show that these long-distance pieces are all correctly taken account of by 
the factor of the parton density in Eq. (^) . 

The scale ^ is the renormalization and factorization scale .|^ Roughly speaking, // sets the 
separation between the parts of the process that we attribute to long- and short-distance 
phenomena. The predictive power of the factorization theorem Eq. @ arises when we set /j, 
to a value of the order of a large physical scale in the problem, say ^ .b,q = ^fac = /^ren) ~ 
IJ-phy ~ \J Then the hard scattering (or Wilson) coefficient a) may usefully be expanded 
in powers of the small coupling as(^). The /U dependence of the parton densities is given 
by the Altarelli-Parisi equation, whose kernel is also perturbatively calculable in powers of 

As explained in itHJ, in the presence of non-zero masses, it is the helicity amplitudes 
which provide the simplest connection between the physical (scalar) structure functions 
Wbn and the corresponding parton-level quantities UJsa- The factorization formula then 
reads: 

WUi(i\^-V) = m®^\a = Y. /^(^,/^)SL(f,-,^,«s(M)), (5) 

where 

= 6W*(g,p)-T^^--6W(g,p) ; = e(,^>{q,k) ■ Cj'^'^ ■ eW[q,k), 

and ei^\q,r) is the polarization vector of the vector boson with momentum q and helicity 
A (= -|-,0, — ) defined with respect to the reference vector r. (Since r is different in the 
definition of W'^ and tD^, the simple relation Eq. (^) is not an obvious consequence of Eq. 
(Q); it follows only because the two reference momenta p and k are collinear, cf. I.) Fig. ^ 
depicts the factorization formula, Eq. (|5[), in a familiar form. 

2.2 Masses and factorization schemes 

The conventional method of calculation of the short-distance coefficients in this and other 
QCD processes is to set to zero the masses of internal lines and external partons in graphs for 

^For simplicity, we do not distinguish between the factorization scale fifac and the renormalization scale 
fJ^ran', wc sct them equal to the same value fi. 
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the partonic subprocesses. Then the resulting infra-red poles (in dimensional regularization) 
are subtracted according the MS scheme. It can be shown that this implies that the ultra- 
violet divergences in the definitions of the parton densities are renormalized by the MS 
scheme also. 

Setting quark masses to zero gives the leading term in an expansion of the short-distance 
coefficient in powers of mg/iiphy This is obviously inapplicable for a heavy quark if we 
want to treat the region where fiphy is not much greater then niQ. However, it is perfectly 
sensible to leave the heavy quark mass in the calculation of Q. We will later show how this 
works in a calculation, and we will verify that the mg — > limit of the coefficient agrees 
with the standard zero mass calculation. The parton densities, including the one for the 
heavy quark, will continue to be defined by the MS scheme. In the case that the heavy 
quark is the charm quark, we will call this the "4-flavor scheme" 

On the other hand, when fiphy *C rn-g, one finds that there are large logarithms of the 
heavy quark mass in all perturbative calculations. That is, the 4-flavor scheme does not 
manifestly exhibit decoupling of the heavy quark. One obvious possibility is to use off-shell 
momentum-space subtractions (which exhibit explicit decoupling) instead of MS. But this 
makes for much more complicated calculations, especially because of the complicated off- 
shell structure of the renormalization counterterms for gauge-invariant operators (such as 
are used to define the parton densities). The method of CWzl§l offers a natural and simple 
way: switch to a "3-flavor scheme" in this region. 

Technically, the CWZ scheme is a hybrid of MS for the light partons and zero-momentum 
subtraction for graphs with a heavy quark line. The scheme has the following advantages: 
It satisfies manifest decoupling, and preserves gauge invariance. The evolution equations for 
the coupling and the parton densities are the same as for QCD with 3 flavors of quark and 
pure MS subtractions. Calculations are quite simple compared with the off-shell scheme. 
Finally, the charm quark density is zero to the leading power in A/rric, so that the charm 
quark mass only appears in Wilson coefficients. 

The 3-flavor scheme is appropriate when fiphy is comparable to or less than about rric- 
When other heavy quarks are present, one defines a series of schemes: 3-flavor, 4-flavor, 
5-flavor etc. The A^-flavor scheme is defined to treat the first flavors of quark as light, 
and the remainder as heavy. It is appropriate when the physical scale of the process, fiphy 
is above the mass of quark and below that of quark -|- 1. We will call N the number 
of active quarks. 

The relation between schemes with different numbers of active flavors is just a case 
of a transformation between different renormalization and faf. tori zation schemes: and the 



matching conditions between the schemes have been calculated ^§1' e3' E3. At the one-loop 
level in the MS scheme, these areE^I just that the coupling and parton densities are the 
same in the two schemes at = mg . Thus a convenient way of implementing them is to use 
3-flavor evolution below fi = rric, to use 4-flavor evolution above that point, with continuity 
at the break point. In this scheme, the use of = mg rather than, say, n = 2mg, is a 
matter of explicit calculation using MS subtraction, and is not a matter of arbitrary choice. 
There are higher order corrections to the matching conditions. Two-loop matching has been 

®Bi_keeping mg non-zero in the Wilson coefficient, the "theoretical inconsistency" described by Gliick 
et al.'B3 does not occur in our approach. See further discussions in Sec. 4.1 



8 



calculated for the coupling, but not yet for the parton densities. 

It is worth noting that existing NLO calculations of heavy quark productionlcJ' E2l 
essentially use the 3-flavor scheme as described above - for all energies, irrespective of the 
order of magnitude of /ip/iy. 



2.3 Contributing Partons and Parton Distributions 

We use the term "variable flavor number scheme" to denote the scheme just described. It is 
implemented^^ by using MS evolution with a number of active flavors that changes as one 
crosses the boundaries \x = mg, where mq is the mass of a heavy quark (charm, bottom, 
etc.). The MT0 and CTEqS parton densities are defined using this method. Thus for a 
given scale /u for the parton densities, all quarks with mass less than [i are treated as partons 
(and thus have associated QCD-evolved parton distributions). For a quark Q with non-zero 
mass rriQ (» Kqcd)-, f^{^i A*) vanishes when \x < mq (i.e., all the heavy quark dynamics in 
this region is in the Wilson coefficients). But when > uiq, /^(^, /i) satisfies the usual MS 
QCD evolution equation (with massless kernel functions) above threshold. Thus, there is 
no fixed restriction on the sum over parton flavor label a in the basic factorization formula 
Eq. (H): depending on the value of the relevant /x (~ ^iphy) of the physical process, the 
correct number of quark flavors appropriate for that energy scale will contribute. 

This conceptual and calculational simplicity has an associated price. In the region just 
above the quark mass (/x ~ thq), defining a parton distribution function for Q with massless 
evolution kernel appears to be somewhat artificial. Indeed the use of f^{^,,lJ-) in a lowest 
order parton model formula for a cross section tends to overestimate the cross section, 
because the parton density does not contain the physical threshold behavior. The errors 
are compensated when one brings in higher order terms in the Wilson coefficient, as we will 
see. Although both schemes are equally correct, in principal, it would seem better to use 
a fixed parton flavor number around threshold, e.g., 3 flavors for charm production. But 
as one goes higher in scale, one is genuinely in the overlap region, where the 3-flavor and 
4-flavor schemes are equally valid. Eventually, the 4-flavor scheme becomes the one which 
describes the underlying physics more accurately. 



2.4 Parton Structure Functions and Hard Scattering Mechanisms 

In the variable flavor number scheme, for both charged current and neutral current pro- 
duction of heavy quarks, initial state quark-partons contribute through the vector-boson 
quark scattering (flavor excitation) subprocess. Fig. ||a and its higher-order corrections; 
whereas the gluon-parton contributes through the vector-boson gluon fusion (flavor cre- 
ation) subprocess, Fig. ^3, and its higher-order corrections. The order quark scattering 
hard amplitude w^g^ is easy to calculate. Since it is obtained from a simple tree diagram, 

we can identify it with the corresponding hard amplitude O'^q which enters the factor- 
ization formula, Eq. (^). In our framework, the explicit results are given in I (Sec. 5 and 
Appendix C). The order a] parton level gluon fusion amplitudes ^j^g'^ are also relatively 
straightforward to evaluate since they are free from singularities when all the quark masses 
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are kept finite.]] 

Note the distinction between the notations cD and co for the structure functions for the 
partonic subprocesses. The unadorned io will be a partonic structure functior^ without 
subtractions, but with non-zero quark masses. In the zero mass hmit uj wih be divergent. 
The hatted quantity lD will have subtractions to remove the infra-red dependence. It is the 
subtracted partonic structure function u) that is to be used in the factorization theorem 
Eq. (I). 

We will present the detailed formulas for lo and u), in the helicity formalism we use, later 
in Sec. ^. Here we focus on the relation between the unsubtracted uj and the subtracted uj at 
order as in order to elucidate the underlying principles. As is well known, this relationship 
is established by applying the factorization formula at the parton amplitude level. This 
provides the exact relation, to this order: 

4'," =E„/,°»s^,„=E.(/;'°'»4':'+/f' (6) 

where we made use of the fact that f^^'^\c) = ^b^{^ ~ 0^ '^Ba'' ~ "^Ba^' ^9 scattering 
begins at order 1. The order quark distribution inside an on-shell gluon, f^^^\ is given 



by 



^5 " 27r ^""m^^^Q' 



where we have used the MS prescription to renormalize the ultra-violet divergences in the 
quark density. Since we have kept a non-zero quark mass, there is no infra-red divergence. 
In this formula, PgQ is the familiar g — > QQ splitting function PgqiO = |(1 — 2^ + 2^'^). 
Eq. (|^) follows from the Feynman rules for parton densities, and it is quite accidental that 
there is no constant term, but only the logarithm times the splitting function. By inverting 
Eq. @, we obtain the formula for the requisite hard amplitude: 

m _,. A(l) _ fQ(l) A(0) 



UJn„ = 



The second term on the right-hand side (henceforth referred to as the subtraction term), 
represents that part of the gluon fusion term which at large ^J.pf^y has the internal quark line 
relatively close to the mass-shell and almost collinear to the gluon. 

2.5 Complete NLO (order a^) Hadron Structure Functions 

Combining these results, we obtain the formula for the physical helicity structure functions 
for heavy quark production on a hadron target: 

W^^ = - E ^ /.'''^"^ ® 42! + ® 41;^ + 0{as'). (9) 



^In the convenliqnal calculation of photo- and lepto-production of heavy quarks in the fixed parton flavor 



number schemec3 these are sometimes called Born terms because they represent the leading order 
contribution to the flavor creation mechanisms. Since the flavor excitation mechanism actually come in with 
one less power of as, we shall state the explicit powers of as to avoid confusion. 

®We make a clear distinction between the concepts of 'structure function' and 'parton density', contrary 
to common usage. 
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Figure 5: Graphical representation of the three terms which enter the master equation, Eq. 
@, for the physical structure functions: the subtraction term is placed in the middle to 
emphasize its similarity both to the quark-scattering (left) and to the gluon-fusion (right) 
contributions. The x on the internal quark line in the subtraction term indicates it is close 
to mass-shell and coUinear to the gluon and the hadron momenta. 

Note that, in the case of a neutral current process, a sum over the two heavy quarks (quark 
and antiquark) in the final state is needed. This is the basic QCD equation for lepto- 
production of heavy quark production in our approach. The subtraction term is placed 
in the middle to emphasize its similarity both to the quark-scattering (left) and to the 
gluon-fusion (right) contributions. On one hand, this term overlaps with the first (quark 
scattering) one due to the common factor uj^q and the approximate equality 

/^(e, /x) ^ f^^^^ ® fff when fi>mQ and a, log ^ ~ 0(1), (10) 

where f^^^^ is given by Eq. I^.f\ On the other hand, its close connection to the last 
(gluon fusion) term originates from Eq. (^). Since it represents the part of the gluon-fusion 
term which is already included in a fully QCD-evolved quark-scattering term, a consistent 
formalism must lead to its subtraction to avoid double counting, as naturally happens here. 
Fig. ^ illustrates the same point graphically. (For clarity, we only show one-half of the cut- 
diagrams for this process, c/.. Fig. ^) The x on the internal quark line in the subtraction 
term denotes the following operation: In the hard scattering part of the middle graph, that 
is, the upper part of the graph, the incoming quark's momentum is replaced by an on-shell 
value with zero transverse momentum. This replacement gives a good approximation when 
the quark is collinear to the gluon, and results in a factor of the order Og distribution of a 
quark in a gluon, Eq. (0). 

The physics behind this formula should be well-known to students of the conventional 
QCD parton model for light quarks.^ However, this formalism has not been invoked in 
existing calculations of heavy quark production. o c3 c3 Rather, they typically use the 
scheme in which the heavy quark only appears in the Wilson coefficients {e.g., a 3-flavor 

^It is straightforward to demonstrate that this expression satisfies the leading order QCD evolution 
equation with the correct boundary condition /^(C, ^^ = nig) = 0. 

^"Although the nature of the subtraction term is usually not transparent to most non-experts since, for 
zero-mass quarks, it is usually identified only as the coefficient of a 1/e pole in the most commonly used MS 
calculational scheme. 
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scheme for charm production and a 4- flavor scheme for bottom production). In that case, 
the quark-scattering (flavor-excitation) contribution is excluded whenever tuq ^ 0, no sub- 
traction from the gluon-fusion contribution is applied in those calculations. There are 
calculations that allow the heavy quark to have a parton density, and thus work well above 
the quark threshold. But these normally set the quark mass to zero in the Wilson coeffi- 
cients and thus are not good approximations when the process is not sufficiently far above 
the quark threshold. {Cf. Ref. g].) 

The variable parton-flavor-number scheme for calculating heavy quark production thus 
represents a natural and correct extension of the usual zero-quark-mass QCD parton frame- 
work to the case of non-zero quark mass. (We avoid using the word "heavy", at least tem- 
porarily, since at this point "light" and "heavy" are relative with respect to the typical 
energy scale in this approach.) This scheme contains all the ingredients of a consistent 
QCD theory of heavy quark production over a wide range of energy scales as mentioned in 
Sec. |l]. In particular, if the initial-state quark (labeled by Q in Eq. (^)) is massive and the 
typical energy scale fXphy is of the same order as mq , then Eq. (lTo|) implies an approximate 
cancellation of the first two terms on the right-hand side of Eq. (P) — thus we arrive at 
W^j^ ~ (g) io^g^ — i.e., dominance of the gluon-fusion mechanism — which reproduces 
the usual picture of heavy quark production in the fixed parton flavor number scheme. (This 
is the region labeled "one large scale {tuq ~ tJ-phy)" in Fig. ^.) 

On the other hand, if either Q is a usual light quark (i.e., thq ~ 0) or Q is massive but 
fJ-phy ^ ^Qi Eq. ( p^ ) does not hold; instead, the subtraction term becomes the dominant 
piece of the gluon fusion contribution (because it has the large logarithmic factor "^^-^ In 

embodying a "collinear divergence"), hence the last two terms almost cancel (leaving only 
a correction term of order with no large logarithm factor) and we obtain W^j^ = 

(g) co^Q + 0{as) - which reproduces the leading order QCD parton model picture 
appropriate for energies much higher than all masses. (This is the region labeled "two large 
scales (1 > mq/iiphy)'" in Fig. ||.) 

Eq- (^) provides a smooth interpolation between the two kinematic regions described 
above, and contains both as special cases. Because the subtraction term represents precisely 
the overlap of the other two, a change in the factorization scale amounts explicitly to a 
reshuffling between the three terms on the right hand side of Eq. (|9|). The differences 
arising from a change in the factorization scale are genuinely of higher order in a^, hence 
are smaller by a factor than if one or more terms are left out.|^ We will demonstrate 
this point in detail in Sec. ^. The differences can be made smaller by using higher order 
terms in perturbation theory. 

Strictly speaking, Eq. ^ is incomplete: we should also add order quark-scattering 

contributions of the form (g) w^g^ — f^^fg ® ^i^Q' • "^^^ ideas are exactly the same 
as discussed above for the order as gluon contributions; but these terms are numerically less 
important because we have many more gluons inside the hadron then sea quarks. (In this 
sense, the order Og Bq scattering terms are effectively one order higher since, for sea quarks, 
~ /at ® is of order as compared to the gluon distribution ffj.) We should also 

mention that Eq. (^) can be generalized to higher orders by the systematic application of 

^^But remember that the gluon distribution is often much larger numerically than a quark distribution. 
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the above scheme. The order hard amphtudes wih be given by formulas generahzed from 
Eq. (P). The differences between these and the corresponding ones already calculated in 
the conventional fixed-number-of-flavor scheme are finite pieces attributable to the change 
of renormalization scheme. 



3 One-loop Gluon-Initiated Parton Structure Functions 



The one-loop forward hard amplitudes for the 2 — > 2 vector-boson-gluon scattering process 



B{q)+9ik) ^qM+Q2iP2) 



(11) 



are given by the cut diagrams shown in Fig. As indicated, we use q to denote the 
momentum of the vector boson, k the momentum of the gluon in the initial state, and 
(pi, P2) the momenta of the quarks in the final state. For flavor changing charged current 
processes (vector boson B = W^), the subscript 1(2) will be associated with the light 
(heavy) quark; for neutral current processes {B = jjZ), both are associated with the heavy 
quark. 



r 



f 



J 



J 





Figure 6: Cut diagrams for order a] vector-boson gluon scattering. 



Since we keep the masses of the quarks (1711,1712) non-zero, these diagrams do not contain 
any singularities: they are infra-red finite since the gluon only appears in the initial state; 
and the internal virtual quark lines never go on mass-shell due to the finite masses. Thus 
it is safe to perform the calculation in 4-dimensional space-time. The parton tensor (cut) 
amplitudes are represented by: 



u.^-^{q,k,p,) =^T.Jdr 



(12) 



where i labels the diagrams; dT2 is the two-particle final-state phase space differential; and 
P1P2 denotes the two propagator factors appropriate for diagram i. 

As mentioned earlier, we find it most convenient to work with helicity structure func- 
tions. They are: 



uj'^iQ'^, s, ml, ml) 



1 



I-' 2 



i'^*{q,k).Nr-ei'\q,k), 



(13) 



where ei^\q,k) is the polarization vector of the vector boson with helicity A (= +,0, — ) 
defined with respect to k. For the general case, the helicity amplitudes exhibit a natural 
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symmetry if we express the vector-boson-quark vertex in terms of the two chiral coupUng 
constants: gn^L- The amphtudes (and Nf") are quadratic in the g's, and we can write 



{K = s,x,a) 



(14) 



with the chiral couphng combinations (symmetric, anti-symmetric, and crossed): 

Cs,a= 9R±gl ; C^ = '29R9L- 

It is convenient to express the results in terms of center-of-mass variables of the parton 
process. We have 



(15) 





= ( 


k, 


0, 


0, 


k 


) 




= ( 




0, 


0, 


-k 


) 




= ( 


El, 


psin^ 


0, 


pcosO 


) 


P'2 


= ( 


E2, 


—p sin 


0, 


—p cos 9 


) 



where the initial state energy and momenta are 

k = (s + Q2)/(2VI) 



and the final state variables are 

P 
El 
E2 



Eq = (s-g2)/(2VI) 

= A(s,m2,m?)/(2^/I) 
(s - m| + mf)/(2\/|) 
{s + ml -mf)/(2VI) 



(16) 



with A (a, b, c) = Va^ + b'^ + - 2ab - 26c - 2ca. 

There are two types of logarithmic terms arising from the phase space integration of the 
propagator factors in the t and u channels, respectively: 

(•s - ni^^ + "'i + A(.s, 111^2- »'i))^ 



Lu — 



log 
log 



Ei+p 
Ei-p 
( E2+P 
\E2-p 



= log 
= log 



Amfs 



(s -I- ml — + A(s, m^, mf))^ 
4mp 



(17) 



3.1 General mass case (Flavor changing charge current interaction) 

For general masses {1711,1712}, the independent right-handed helicity structure functions are: 



UJ 



+ — 



+ 



+ 



Lt 
+ Lu 



1 

2 k \ k 

1 Eo ( Eo 

- + — — - 1 

2 k \k 



{Q^ + 1)2 

2 k \ k 2sfc2 



1 £2 

2 k [ k 



1 + 



m^mf — ml) 

2iP 



+ 



2p 

71 



8mim2pV§ 

(g2 + s)2 

p{m2 — mf ) 
k'^s/i 



(18) 
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where the subscripts {s,x,a) refer to the chiral combinations of Eq. (0). 
The longitudinal helicity ones are: 



Wo 



-{Lt + Lu) 



{ml + m|)(Q^ - {ml - m{Y - 2sk^) 



(771-2 ~ rnf) \ m^mf 



+ 



+ 



-{Lt — Lu 
P 



^ Eg{ml - ml) 
k'^\/I 

[{ml - mlf + Q2(2Q2 _ (^2 ^ ^2))] 



-{Lt + Lu) mim2 



1 s — {ml + m^) 
Q2 + 



+ 



2mim2P 
A;2-y| 



0. 



(19) 



And the left-handed helicity structure functions are related to the right-handed ones 
by the symmetry relations: 



UJ^ = UJ' 



-W ' 



(20) 



3.2 Equal mass case (Flavor non-changing neutral current interaction) 

When the masses are equal, mi = m2 = m, we obtain, as a special case of the above: 

\/I -|- — 4r?T,2 



and 



where 



LOX. 



Lt = Lu = L = 2log 



2m 



{ Q^ + s^) _ (s-Q2)2A 

,2 + ^)2 

4m2(s — 2m2'' 



'^^ ^ (Q2 + s)2 s{Q2 + s)2 



Am^/S. 



-L 



(Q2 + 5)2 (Q2 + J)2 

2m2(4m2Q2 + 3g4_4g2^_^2) 4(Q2 _ ^^2-,^ 



g2(Q2 + 1)2 



+ 



(Q2 + s)2 

2m2(-4m2(52 + g4_^4Q2^_^^2) 4^2^ 
-tj „r. , „o ttt; h 



g2(Q2 + ^) 



A = A(s, m ,m ) = y s{s — Am^). 



(Q2 + s)2 



(21) 



4 Mass-singularities, Collinear Divergences, Subtractions and 
Infrared-safe Amplitudes 

In the fixed parton-flavor-number calculational scheme, the results of the last section rep- 
resent the full answer to the vector-boson-gluon fusion production of heavy quarks at the 
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order level. These results contain terms which become large as one (or both) of the 
quark masses are small compared to the characteristic energy scale. These mass singular- 



ity terms are isolated by taking the mi (m2) — > limits of Eq. (18-2C). They arise from 
the configuration in the phase space integration when an internal quark line becomes al- 
most on-shell0 and collinear to the initial state gluon. In our scheme with scale-dependent 
parton- flavor-number, these terms are included in the quark scattering contribution with 
properly evolved quark parton distributions. As explained in Sec. |2.5| the QCD formalism 
provides a natural procedure to subtract these terms from the gluon-fusion amplitudes to 
avoid double counting. In this section, we identify the subtraction terms in some detail. 

4.1 Unequal Mass (flavor-changing) Case (mi 0) 

In the limit mi 0, we drop all mi dependencies in Eq. ([l^pO|) except inside the logarithm 
where the mass singularity resides. We also replace m2 by mq to emphasize that the 
remaining mass is associated with a quark that is heavy in absolute terms (i.e., compared 
to A). The arguments of the t- and u-channel logarithmic factors become 



Ei+p {Ei+pY 4p' 

Lt = log — = log 2 ^ — 

El — p m^ mi 

Lu : log — = log 2 — 

E2 — p m2 mp 



(22) 



where the magnitude of the quark 3-momenta is given by p = {s — mg) /{2s). We note 
that Lt is the only factor which contains the mass singularity associated with mi ^ . It 
arises from the collinear integration region of the t-channel propagator factor over the quark 
transverse momentum, Fig. ^a,b,c. This is seen as follows: in the limit under consideration, 
the singular factor in the integrand of Eq. ( p!3[) is l/(t — mf) oc l/{pt + m^). We can take 
the mi — > limit everywhere except in this factor where it must be retained to cut off the 
collinear (i.e. pt = 0) singularity. The leading behavior is obtained by keeping the constant 
term of the Taylor expansion of the numerator function in p^. We obtain therefore, 

P?max dpf max r , rt f 1 \ ^OQ\ 

2 = log 2~ = + 0{l) (23) 

pf + mf mf 

where we used p = pt ma.x for the quark lines, and explicitly displayed the role of ml as the 
cutoff for the collinear singularity. 

Isolation of Infrared-sensitive terms: 

We can isolate the mass (collinear) singularity from the (process-dependent) dynamics 
by writing 

Lt = log ^ + log^%S£ (24) 

where we have introduced an arbitrary scale parameter fif — the factorization scale — 
which separates the low pt collinear region from the true hard scattering region. For an 
appropriately chosen fj,f {e.g., some external physical scale independent of s), the first term 

^^Relative to Q^. 



16 



on the right-hand side (the colhnear term) contains the mass singularity which is to be 
subtracted from the gluon-fusion contribution (to make it infra-red safe) and resummed 
into the QCD evolved quark distribution function in the quark-scattering contribution. 



The infra-red sensitive terms can be collected by substituting Eq. (|24| ) in Eq. (|18| -| 
We obtain the following non-vanishing amplitudes: 



1 • Pg^g(Xm)l0g ( 



m 



^i""'^ = 1 • Pg^,{x„^)\0g(^ (25) 



m 



m 



2 / ,,2 



with Xm denoting the scaling variable, 



2 




g2 Q2 



2k -q s + Q^ 
and Pg^q the usual gluon to quark splitting function, 



(26) 



P,^,(x) = i[(l-x)2 + x2] (27) 

They are exactly proportional to the leading order quark scattering amplitudes in the same 
limit. Comparing Eq. (p5|) with the order result from paper I, we obtain: 



„2 

u'Jl ^ = J^f (mi = 0) • ^Pg^q{xm)9{i^f - mi) log ^ (for all A) (28) 

This is anticipated in the discussion of the factorization theorem in Sec. |2|, cf. the second 
term on the right-hand side of Eq. (^) . 

Mass (Collinear) Subtraction and Infrared-safe Amplitudes: 

As explained in Sec. ^, the collinear configuration from which the above infrared sensitive 
amplitudes originate corresponds physically to the overlapping region of the quark scattering 
and gluon fusion production mechanisms. We need to subtract these amplitudes from uj^g 
to avoid double counting and to obtain infra-red safe hard gluon scattering amplitudes 
{cf. Eq. (^)). We have a certain freedom in choosing the subtraction terms: as long 
as they contain the leading infrared sensitive terms identified above, any specific choice 
made for these overlapping terms effectively defines the factorization scheme inherent in 
the master (factorization) equation, Eq. (^). However, this freedom must be constrained 

by two consistency requirements. On the one hand, the order quark amplitude uJj^^^ 
should be identical in its two occurrences in Eq. (^). On the other hand, the logarithmic 
term in Eq. ( p8|) must agree with the scheme used to define parton densities. 
There is a certain amount of choice possible. 

We consider the case of unequal masses, one large and one small: m2 /U/ = fiphy ^ 
mi » Aqcd- Then the parton density that is needed for quark 1 is defined by MS. Consider 
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^i^qi ' "^hich is given by the Born graph for scattering off the relatively light quark 1 to make 

the heavy quark 2. There are two obvious choices for the value of mi in oj^'^ : either replace 
the mass mi of the lighter quark to zero, or leave it at the physical value. When mi is small 
compared with the physical scale for the hard scattering, ^phy, it certainly does not matter 
what we do, since Wilson coefficients have a finite zero-mass limit. But our formalism also 
extends to the region where mi is not negligible. Then the physics is correctly given by the 
order ag graph for Bg fusion, the third term i^^^^^ in Eq. @ — with the mass mi kept at 
its physical value. (This is the region where the 3-flavor scheme applies and it treats mi as 
"heavy".) For consistency, it is then advantageous to keep mi also non-zero in the order 
Wilson coefficient ^ even if this is not absolutely required by the formalism — the first 
two terms in the master equation are guaranteed to cancel near threshold (c/. Sec. 2^) as 
long as the same choice of is made in both terms. 
Thus we shall define the subtraction term as 



, ,Sub A 
^Bg 



"T'l)log;^ 



(for all A). (29) 



with the full mass dependencies in uj^^\s,mi,mQ), and with the understanding that the 

same choice of oj^^^ is to be made in the order quark-scattering term in the master 
equation, Eq. (|9|). The choice of the MS scheme for the parton densities implies the precise 
formula given in Eq. ( p9[) , which has a logarithmic term, but no constant term. This follows 
from a calculation of the one-loop density of a massive quark in an on-shell gluon. 

The properly subtracted, infra-red safe hard gluon scattering amplitudes are then given 

by 



ujBg{s,Q ,mi,mQ,fif) = uJBgis,Q ,mi,mc 



ujl^f^{Q^,mi,mQ,Hf) 



(30) 

with uj^g{s,Q'^,mi,mQ) given by Eqs. (|l8|)-(|0l) of Sec. |T| and uj^^f ^{Q"^ ,mi,mQ, ^f) by 



Eq. (|29D above respectively. 

The explicit expressions for tDg^ in the general case is not particularly illuminating. We 
give below the results in the mi — > limit: 



Pa 



log 



4p: 



t max 



x{l-2xf {s-ml) 



/x2 



(1 



9~*q 



. Xr, 



+ 2x{l - 2x) 



Pg^g{Xm)'iOg 



4p- 
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(2x2) 



m} 



Pg^g{Xm)+2x{l 
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2i,' 




(31) 
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llog^%2£ 
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m 



+ x^ 
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Of course, for energy scales much larger then m2 (w-q), there is also a mass singularity 
associated with the factor log(m2//Uj) which resides in Lu- In our approach, an infra-red 
sensitive term analogous to Eq. ( p8| ) with mi <-> 1712 should then be subtracted from the order 
Os gluon-fusion amplitudes. The subtracted part, again, is included in the corresponding 
order quark scattering amplitude — with incoming quark "2" — which represents the 
resummed result of all such terms to arbitrary orders. The resummation is performed by 
the Altarelli-Parisi evolution of the parton densities. 



4.2 Equal Mass Case (mi = m2 = m — 0) and Comparison with MS 
scheme results 

Corresponding results for the case of flavor non-changing neutral current interactions can 
be obtained from the above general results by setting mi = m2 = m and choosing the 
appropriate couplings. We give a few explicit formulas for illustrative purpose and for 
establishing the relation of our subtraction scheme to the MS scheme. For the equal mass 
case, the infra-red sensitive logarithm factor is 

L = log = log 4 = l°g 4 + log 4 

m^ m^ m^ fij 

In the m ^ limit, the non- vanishing helicity amplitudes are (keeping m only in the 
other-wise divergent logarithm) 



u;J = 2Pg^g{x) log — - (1 - 2x) 
m^ 

LOg = 4(1 - x)x 



(32) 



After subtracting the mass-singularity, we obtain the infra-red safe hard amplitudes for zero 
mass quarks 

2+ = Q~ = 2Pg_,,{x) log 4 - (1 - 

= 2P,_,(x) (log^ + logi^)-(l-2x)2 (33) 



cD° = 4(1 - x)x 



X 



These correspond to the "Wilson coefficients" for deep inelastic scattering, as usually cal- 
culated in the MS scheme. It is straightforward to verify that they, indeed, are identical to 
the MS Wilson coefficients. IHl Hence, our subtraction prescription (applicable to the gen- 
eral mass case) reduces to the MS-scheme of subtraction of collinear singularities in the 
zero-mass limit. This is simply a consequence of our choice of the MS scheme to define the 
parton densities, which then resulted in Eq. (0) for the one-loop value of the quark density 
in an on-shell gluon. A change of definition of the parton densities would have added an 
infra-red safe term to Eq. (0), and there would be corresponding terms to be added to the 
other terms in the formulae we have written. 
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5 Results on Structure Functions 



We shall now study the numerical significance of the quark scattering (QS) mechanism 
compared to gluon fusion (GF) in this unified framework. For the case of charged current 
production of a heavy quark from a light quark, we have demonstrated in a previous pub- 
lication that the two basic processes are of the same size numerically; hence a quantitatiw 
QCD analysis must incorporate both in a consistent manner such as formulated above. 1^ 
In the following, we shall concentrate on heavy quark production by neutral current in- 
teraction which is of great interest both at fixed target and ep collider energies. For the 
purpose of this paper, we shall present results on the heavy quark production structure 
functions F2{x,Q) (which can be applied at all relevant energies). Phenomenological re- 
sults on cross-sections for various interesting processes at specific energies of fixed target 
and HERA experiments will be pursued in a subsequent study. For practical reason, we 
show mostly results on charm production. At the end, we also show some corresponding 
results on b-production. 

Two sets of parton distributions are used in the following study: the "next-to-leading 
order (NLO) set" consists of the CTEQ2M distributions, H and the "leading order (LO) 
set" which is generated from initial distributions at Qq = 1.6 GeV taken from CTEQ2M 
but evolved using LO evolution kernel only. The differences between using LO and NLO 
are quite substantial, as regards the charm-quark distribution, and this is a symptom that 
higher order terms are important. 



5.1 Scale-dependence of the Structure Function 

As with all applications of the perturbative QCD parton formalism, a theoretical uncer- 
tainty about these calculations concerns the choice of factorization scale (which we identify 
with the renormalization scale, c/., Sec. |lj). Whereas some scale dependence of the theo- 
retical prediction is unavoidable to any given order in Og, an excessive sensitivity to the 
scale parameter usually signals a large theoretical uncertainty. This has been a worry 
for existin, 
scattering. 



) calculations of charm and bottom production, especially for hadronic 
I In order to make clear how the physical results presented later might 
depend on the (implicit) choice of scale, we first investigate the scale dependence of the var- 
ious contributions which enter our calculations. In particular, we demonstrate that whereas 
both the quark-scattering (QS) and the gluon-fusion (GF) terms show substantial scale 
dependence, these dependences are opposite in direction and they compensate each other 
when combined according to the variable flavor number scheme, cf. Eq. (|9|). 

Fig. 0a shows F2{x, Q; /x) as a function of // for x = 0.05 and Q = 10 GeV, calculated 
using the LO parton distribution set. We display the //-dependence of the QS (long-dashed 
line), the GF (dotted line), and the Subtraction (short-dashed line) terms individually 
along with the combined (solid line) Total result. We see that the rapid rise of the QS 
and the Subtraction terms together with the somewhat gentle fall of the GF contribution 
combine to make the Total result substantially more stable then either of the two individual 
production mechanisms. We also note the following important features of Fig. for n 



below the heavy quark mass, the QS and subtraction terms vanish by definition {cf. Sec. 2^ 



and Eq. (|T^)), we have -F^J""* = -^2^^! ^'^^ J^^t above the mass threshold, the QS and 
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Figure 7: Scale dependence of the contributing terms to F2{x,Q). The factorization scale 
fi is shown in units of the physical scale 2Mq. 



Subtraction contributions nearly cancel according to Eq. (10), we have -^2"°* — -^2'^ (and 
the difference has reduced fi dependence); but for fi much larger than the mass threshold, QS 
(after subtraction) makes a substantial contribution to -FJ""*, and -F^'^ ceases to be a good 
approximation. These features provide supporting evidence to the theoretical discussion of 
Sec. |2.1[ - p.5| , as will be reinforced by results to follow. 

Fig. ^b shows the corresponding results obtained with the NLO parton distribution set. 
All the qualitative features are the same as above. However, the cancellation between the 
QS and Subtraction terms above the heavy quark threshold is not as complete as in Fig. ^a. 
(The reason is that Eq. (10) is not as good an approximation as in the previous case, because 
the parton densities are evolved at the NLO, but Eq. (10) is only used at order a] in our 
calculation.) Hence the contribution of the QS production mechanism (after subtraction) to 
the structure function is even more significant for all values of ^ above the threshold. The 
NLO parton distribution functions contain resumed sub-leading logarithms of fi/niQ, thus 
the difference between the QS and Subtraction terms contains additional pieces of higher 
order terms not present in the order a] GF calculation included in this study. Note the 
complete stability of the Total curve against the choice of fj, in this case. 

These results imply that, in subsequent discussions of the x- and Q-dependencies of 
F2{x, Q), the choice of the scale /i can shift the QS, GF and Subtraction terms individually 
by considerable amount, but it will not affect the Total answer by nearly as much. This fact 
underscores the intrinsic inter-dependency of the two heavy quark production mechanisms 
(as a basic quantum mechanical mixing effect). It also ensures that we actually have a fair 
range of freedom of choice of the scale — which we can take advantage of whenever there 
is good physics reason to do so. 



5.2 Choice of Scale 

When Q ^ rnq, the natural hard scale of the production process is of the order Q. For 
Q ~ 0{mQ), jjL can in principle be any combination of Q and mq which is of the same 
order of magnitude. To make an intelligent choice however, it is important to be guided 
by relevant physical considerations. Since we know that gluon fusion represents the correct 
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physics near the threshold for heavy quark production, it is appropriate to choose a /i such 
that the quark scattering contribution (along with the subtraction term) becomes small 
in this region. It also makes sense to let the latter vanish when Q < mq. As a concrete 
example, the following ansatz for the scale /i satisfies all these requirements. 

^2 =rnl + cQ^l-^r for Q>mQ (34) 

= rriQ for Q < mq 

The results presented below are obtained using this ansatz with c = 0.5 and n = 2. Of 
course, as is the case with all pQCD calculations, an infinite number of other choices are 
also acceptable. However, the results of the last subsection ensure that, in our formalism, 
qualitative features of the answer will be common to most reasonable choices — as we have 
verified by actual calculation with a variety of prescriptions for /x. (This is obviously not 
the case if the QS or GF production mechanisms are taken individually, particularly if an 
extended range of Q is involved.) 

5.3 Behavior of Structure Functions and the Interplay between the QS 
and GF Production Mechanisms 

Fig. ^a shows F2{x, Q) for charm production as a function of Q for fixed x = 0.01 using LO 
parton distributions. The lines are labeled the same way as in Fig. |^ For our particular 
choice of scale, the QS contribution (long-dashed line) emerges from threshold and becomes 
comparable in size to the GF contribution (dotted line) beyond around 5 GeV. Because the 
subtraction term tracks the QS term rather closely throughout the kinematic range except 
at the very large Q end, the net contribution of these two is quite negligible, hence the 
Total curve (solid line) stays very close to the GF one except for very large Q. 




10' 10' 
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Figure 8: Q dependence of F2{x,Q) at a; = 0.01 and x = 0.1 calculated using LO parton 
densities. 

Fig. |p3 is analogous to Fig. ^a except that we now look at the behavior of the various 
contributions to F2{x,Q) at a larger x value, x = 0.1. The picture is somewhat different 
and, for the purpose of illustrating the physics underlying our approach, more illuminating. 
For the same choice of scale as above, the QS contribution rises rather steeply, and overtakes 
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the GF contribution almost immediately above threshold. Of special interest is the behavior 
of the subtraction term (short-dashed line): it tracks the QS contribution above threshold, 
as noted before, then turns to follow the GF curve at large Q. The latter behavior follows 
from the definition of the subtraction as the leading collinear log term of the (order Us) GF 



contribution, cf. Sec. 2.5 and Sec. 4.1. As a result, the complete QCD result (the "Total" 
curve) follows the GF term at low Q, but it approaches the QS contribution at high Q. 
We note, this same behavior is also present in Fig. |8|a, for x = 0.01. However, there it is 
obscured by the closeness of the QS and GF terms due to our particular choice of scale and 
parton distributions. 
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Figure 9: Q dependence of F2{x,Q) at x = 0.01 and x = 0.1 calculated using NLO parton 
densities. 

In Figs, ^a and §3, we show the same curves as in Figs, ^a and |8|b respectively, now 
calculated using NLO parton distributions. The same physics effects are clearly displayed 
in both cases, only in more dramatic proportions. 

The general feature of close interplay between QS and GF production mechanisms follows 
directly from our basic premises discussed in the introductory sections. Sec. || and Sec. |^: 
when the relevant physical energy scale (Q) is comparable to the mass of the heavy quark 
(mg), this quark behaves more like a heavy particle rather than a parton, hence GF is 
the dominant production mechanism but when Q ^ mg, it behaves characteristically 
like a light quark parton almost by definition, and QS becomes the dominant process. 
This intuitively reasonable behavior naturally emerges from the variable (scale-dependent) 
flavor-number scheme of calculating massive quark production, Eq. (^). 

We now look at F2{x, Q) as a function of x at fixed Q. Figs, and b show the various 
terms for Q = 10 GeV calculated using LO and NLO parton distributions respectively. 
We see the characteristic rise of the structure function toward small x. For fixed Q, the 
large x limit corresponds to the total final-state energy W Wthreshoid (for heavy quark 
production); all contributions become small. As x — > 0, one moves away from the threshold 
region, the relative size of the various terms are sensitive to the choice of scale and the 
choice of parton distribution functions. 

^^This is the basic tenet of the analysis of Collins, Soper and Sterman, Ref. The underlying physics 
was quantitatively demonstrated previously in a different process by Olness and Tung in Ref. p9[|. 
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These results clearly illustrate the importance of the QS ("flavor excitation" in old 
literature) mechanism for heavy flavor production. The GF ("flavor creation") mechanism 
provides a natural explanation of the production of the heavy quarks not far above the 
threshold; but it is not adequate to account for this process when the energy scale becomes 
large. Part of the QS contribution (with proper subtraction) will be included when the next 
order (i.e., a^) GF hard-scattering terms are included. |p6[ However, the latter only contain 
terms to order In ^(^/Mq), whereas the QS contribution represents the resummed results 
of all such terms to arbitrary orders. 



5.4 B-quark Production 

Results on b-quark production are similar to those shown for charm. We show only one 



plot. Fig. 11, which shows F2{x,Q) as a function of Q for fixed x = 0.01. We see that the 
features are entirely similar to those seen in Figs. |8| and ^ for charm production. 



0.12 




b production 




0.1 


Total 

■ QrkSc 






0.08 


■ Subtr 

GluFu 






0.06 








0.04 






x = 0.05 


0.02 






LO PDF 














10' 


Q (GeV) 







b production 






Total 






QrkSc 






Subtr , •- 






GluFu ' / 








X = 0.05 






NLO PDF 









Q (GeV) 



Figure 11: Structure function F2(x,Q) for b-production vs. Q at x = 0.01. 
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6 Discussions 



We have shown in this paper that the currently available fixed-order QCD calculations of 
heavy quark production using the flavor creation mechanism alone have to be generalized to 
include flavor excitation with scale-dependent number of quark-flavors in order to account 
for the appropriate underlying physics at all energies. A consistent scheme to implement 
this generalization is formulated in detail to order Ug for leptoproduction here. The method 
can be extended to higher orders: in addition to the flavor-creation diagrams which are 
already calculated in the literature, one needs to add order ag flavor-excitation (vector- 
boson scattering off "heavy quark" partons) contributions and perform the appropriate 
subtraction. As remarked earlier, these two contributions are numerically comparable in 
spite of the formal difference in the power of as by one. Exactly the same principles apply 
to hadroproduction. As the results of the previous section show, the inclusion of the right 
physics in these calculations can be expected to improve the theoretical accuracy of QCD 
predictions — as illustrated by the reduced dependence on the choice of (spurious) scale. 

The existing fixed order calculations have a natural region of validity: when ^ip^y is 
of the same order of magnitude as mq and x is not too small. Our proposed scheme for 
heavy quark production calculation contains the right physics when fiphy becomes much 
larger than mq. An important question to ask is then: where exactly lies the transition 
region and what formalism should be used in the transition region? In order to discuss 
this question in specific terms, let us choose the case of charm production, ignoring b- and 
t-production completely. 

Just above the threshold for producing charm, we all agree that fixed-order calculations 
using flavor creation (gluon-fusion plus light-quark scattering) alone should be reliable — 
provided that the parton distribution functions used in the master formula are generated 
in using evolution kernels with effective flavor number n^ff = 3! We call this method of 
calculation the 3- flavor scheme. (Note, this is not usually done in the published literature. 
Rather, authors of the existing calculations invariably use canned parton distributions con- 
taining 4 or 5 quark-partons, depending on the energy scale. This procedure is, in princi- 
ple, inconsistent.) Far above threshold, say fipf^y > 20m,Q, our proposed scheme (including 
charm as one of the partons) should become the more reliable method. Here, we must use 
parton distributions generated with rieff = 4 QCD evolution equation. We call this the 
4-flavor scheme. The question is: what scheme should be used in the transition region, say 
10 mg > fJ^pfiy > 20 niQ? Although this is one of those elusive questions in pQCD which 
defies definitive answer, our best answer is: the region of transition should be considered 
as the region of co-existence of the two scheme — both schemes should be close to the real 
answer and either one can provide a reasonable result. For this view point to be viable, 
the answers obtained from the two approaches must be close to each other in the transi- 
tion region, with the difference being of the order of the next order of perturbation theory, 
without large logarithms. The example given in the previous section indicates that indeed 
this is the case for some range of value of Uphy. In fact, in the absence of reliable prediction 
of where the transition region lies (as with "when should Bjorkcn scaling sets in?"), the 
requirement of approximate equality of the prediction of the 3-flavor and 4-flavor scheme 
gives the most reasonable criterion for identifying where the transition takes place. (In the 
above discussion, we used 10 mq > fi^f^y > 20 mg only as an illustration. The appropriate 
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numbers should be found in this phenomenological way.) 

The present paper focuses on the motivation and the physics ideas. A detailed compar- 
ison of results from the two schemes and a study of how does the transition takes place, 
as well as phenomenological results pertaining heavy quark production at fixed-target and 
collider experiments will be presented in subsequent publications. 
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